MULTIPLICATIVE FUNCTIONAL FOR REFLECTED BROWNIAN 
MOTION VIA DETERMINISTIC ODE 



KRZYSZTOF BURDZY AND JOHN M. LEE 

Abstract. We prove that a sequence of semi-discrete approximations converges to a multi- 
plicative functional for reflected Brownian motion, which intuitively represents the Lyapunov 
exponent for the corresponding stochastic flow. The method of proof is based on a study of 
the deterministic version of the problem and the excursion theory. 



1. Introduction 

This article is the first part of a project devoted to path properties of a stochastic flow 
of reflected Brownian motions. We will first outline the general direction of the project and 
then we will comment on the results contained in the current article. 

Consider a bounded smooth domain D C R n , n > 2, and for any x € D, let Xf be 
reflected Brownian motion in D, starting from Xq = x. Construct all processes X x so that 
they are driven by the same n-dimensional Brownian motion. It has been proved in [BCJ] 
that in some planar domains, for any x 7^ y, the limit lim^oo log \Xf — Xf\/t = A(D) exists 
a.s. Moreover, an explicit formula has been given for the limit A(D), in terms of geometric 
quantities associated with D. Our ultimate goal is to prove an analogous result for domains 
in R n for n > 3. 

The higher dimensional case is more difficult to study for several reasons. First, we believe 
that the multidimensional quantity analogous to A(D) in the two dimensional case cannot 
be expressed directly in terms of geometric properties of D. Instead, it has to be expressed 
using the stationary distribution for the normalized version of the multiplicative functional 
studied in the present paper. Second, non-commutativity of projections is a more challenging 
technical problem in dimensions n > 3. 

The result of |BCJ] mentioned above contains an implicit assertion about another limit, 
namely, in the space variable for a fixed time. In other words, one can informally infer the 
existence and value of the limit \im. £ ^(Xf +£V — Xf) / e = At"v, for v e R n . The limit operator 
At, regarded as a function of time, is a linear multiplicative functional of reflected Brownian 
motion. Its form is considerably more complex and interesting in dimensions n > 3 than in 
two dimensions. 
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Our overall plan is first to prove the differentiability in the space variable stated in the last 
paragraph. Then we will prove the existence and uniqueness of the stationary distribution for 
the normalized version of A%. And then we will prove the formula for the rate of convergence 
of \X*-X\\ to 0, as t -> oo. 

The immediate goal of the present paper is much more modest than the overall plan 
outlined above. We will deal with some foundational issues related to the application of 
our main method, excursion theory, to the convergence of semi-discrete approximations to 
the multiplicative functional described above. We will briefly review some of the existing 
literature on the subject, so that we can place out own results in an appropriate context. 

The multiplicative functional A t appeared in a number of publications discussing reflected 
Brownian motion, starting with [XJ HKlj . and later in \IK2\ |Hj . None of these publications 
contains the analysis of the deterministic version of the multiplicative functional. This is 
what we are going to do in Section [2j In a sense, we are trying to see whether the approach 
of [LSj could be applied in our case — that approach was to develop a deterministic theory 
that could be applied to stochastic processes path by path. Unfortunately, our result on 
deterministic ODE's do not apply to reflected Brownian motion, roughly speaking, for the 
same reason why the Riemann-Stiltjes does not work for integrals with respect to Brownian 
motion. 

Nevertheless, our deterministic results are not totally disjoint from the second, probabilis- 
tic section. In fact, our basic approach developed in Lemma 12.91 is just what we need in 
Section [3j The main theorem of Section [3] proves existence of the multiplicative functional 
using semi-discrete approximations. The result does not seem to be known in this form, 
although it is obviously close to some theorems in [A] HKl} |H] . However, the main point 
is not to give a new proof to a slightly different version of a known result but to develop 
estimates using excursion techniques that are analogous to those in [BCJ], and that can be 
applied to study Xf — X\. 

We continue with some general review of literature. The differentiability of Xf in the 
initial data was proved in [DZ] for reflected diffusions. The main difference between our 
project and that in |DZj is that that paper was concerned with diffusions in (0, oo) n , and 
our main goal is to study the effect of the curvature of dD. Deterministic transformations 
based on reflection were considered, for example, in \LS\ IDIl [DRj. Synchronous couplings of 
reflected Brownian motions in convex domains were studied in |CLJ1[ ICLJ2] , where it was 
proved that under mild assumptions, Xf — Xf is not at any finite time. Our estimates in 
Section [3] are so robust that they indicate that Theorem 13 . 21 holds for the trace of a degenerate 
diffusion on dD, defined as in [CSl [MO] , with the density of jumps having different scaling 
properties than that for reflected Brownian motion. In other words, the main theorem of 
Section [3] is likely to hold in the case when the trace of the reflected diffusion is any "stable- 
like" process on dD. We do not present this generalization because, as far as we can tell, 
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the multiplicative functional A t does not represent the limit lim e | (Xf +ev — Xf)je for flows 
of degenerate reflected diffusions. 

We are grateful to Elton Hsu for very helpful advice. 

2. Deterministic differential equation 

2.1. Geometric Preliminaries. Throughout this section, M will be a C 2 , properly em- 
bedded, orientable hypersurface (i.e., submanifold of codimension 1) in R n , endowed with a 
unit normal vector field n. The properness condition means that the inclusion map M <— ■> R n 
is a proper map (the inverse image of every compact set is compact), which is equivalent to 
M being a a closed subset of R n . For any R > 0, let Mr denote the intersection of M with 
the closed ball of radius R around the origin in R n , and note that M R is a compact subset 
of M. 

We consider M as a Riemannian manifold with the induced metric. We use the notation 
(•, •) for both the Euclidean inner product on R n and its restriction to T X M for any x G M, 
and | • | for the associated norm. 

For any x G M, let ir x : R n — > T X M denote the orthogonal projection onto the tangent 
space T X M, so 

7r x z = z - (z,n(x))n(x), (2.1) 

and let S(x) : T X M — > T X M denote the shape operator (also known as the Weingarten map), 
which is the symmetric linear endomorphism of T X M associated with the second fundamental 
form. It is characterized by 

S(x)v = -d v n{x), v G T X M, (2.2) 

where d v denotes the ordinary Euclidean directional derivative in the direction of v. If 
7 : [0, T] — > M is a smooth curve in M, a vector field along 7 is a smooth map v : [0, T] — > M 
such that v(t) G T^M for each t. The covariant derivative o/v a/ong' 7 is given by 

X> t v(f) := v'(0 - (v(0,5(7(0)y(0>n(7(*)) 
= v'(t) + (v(t),^(no 7 )(t))n(7(t)). 

The eigenvalues of <S(x) are the principal curvatures of M at x, and its determinant is the 
Gaussian curvature. We extend S(x) to an endomorphism of R™ by defining S(x)n(x) = 0. 
It is easy to check that S(x) and ir x commute, by evaluating separately on n(x) and on 
v G T X M. 

The following lemma expresses some elementary observations that we will use below. 
Most of these follow easily from the fact that smooth maps satisfy uniform local Lipschitz 
estimates, so we leave the proof to the reader. For any linear map A: R™ — > R n , we let ||^4|| 
denote the operator norm. 
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Lemma 2.1. For any R > and T > 0, there exists a constant K depending only on M , 
R, and T such that the following estimates hold for all x, y G Mr, < / < T , t > and 
z G R n : 



IKx - %l 


<K\x-y\. 


(2.3 


l|5(x)|| 


< K. 


(2.4 


l|5(x) -5(y)|| 


<K\x-y\. 


(2.5 


|| e *5(x)|| 


<e Kt . 


(2.6 


le^zl 


>e- Kt \z\. 


(2.7 


||e' 5(a:) - Id || 


< Kl. 


(2.8 


e i5(x) _ e lS(y)u 


<Kl\x-y\. 


(2.9 


\n(x) - n(y)\ 


<K\x-y\. 


(2.10 



Another useful estimate is the following. 

Lemma 2.2. For any R > 0, there exists a constant C depending only on M and R such 
that for all w, x,y,z G Mr, the following operator-norm estimate holds: 

\\ir z o (-Ky - Tr x ) o tt w \\ <C(\w- y\ \y - z\ + \w — x\ \x-z\). 

Proof. Using the fact that n is a unit vector field and expanding \xi{x) — n(y)| 2 in terms of 
inner products, we obtain 

(n(z), n(y)) = |(|n(x)| 2 + |n(y)| 2 - |n(.x) - n(y)\ 2 ) = 1 - ||n(x) - n(y)\ 2 . 

Suppose w, x,y, z G Mr and v G R n . If n w v = 0, then the estimate holds trivially, so we 
may as well assume that v G T W M. Expanding the projections as in (12.11) and using the fact 
that tt w v = v, we obtain 

= tt z (v - (v, n(y))n(y)) - tt z (v - (v, n(x))n(x)) 
= (v - (v,n(y))n(y)) 

- ((v,n(z)}n(z) ~ (v,n(y))(n(y),n(z))n(z)) 

— (v — (v, n(x))n(x)) 

+ (( v > n ( 2; )) n ( 2 ) - (v, n(x))(n(x),n(2;))n(z)) 
= -(v,n(y))n(y) + (v,n(x))n(x) 
+ (v,n(y))(l-±\n(y)-n(z)\ 2 )n(z) 

— (v, n(x)}(l — ||n(x) — n(z)\ 2 ^n(z) 

= -(v, n(y))(n(y) - n(z)) + (v, n(x))(n(x) - n(z)) 

- l( v > n (2/))l n (2/) -n(z)\ 2 n(z) + |(v, n(x))|n(x) -n(z)| 2 n(z). 



Multiplicative functional for reflected Brownian motion 5 

Using the fact that (v, n(w)) = 0, this can be written 

ix z (ix y - 7rx)-K w v = -(v, n(w) - n(y))(n(y) - n(z)) 
+ (v, n(w) — n(x))(n(x) — n(z)) 
- |(v,n(w) - n(y))\n(y) - n(z)\ 2 n(z) 
+ |(v, n(w) — n(x))\n(x) — n(z)\ 2 n(z). 
The desired estimate follows from (12.101) and the fact that 

\n(x) — n(y)\ 2 < (\n(x)\ + |n(y)|) \n(x) — n(y)\ < 2K\x — y\. 

□ 

2.2. Analytic Preliminaries. Let T be a positive real number. We let BV([0, T]; R) denote 
the set of functions u : [0, T] — > R of bounded variation, and NBV([0, T}; R) C BV([0, T}; R) 
the subset consisting of functions that are right-continuous. By convention, we will consider 
each u G NBV([0,T]; R) to be a function defined on all of R by setting u(t) = for t < 
and u{t) = u(T) for t > T; the extended function is still right-continuous and of bounded 
variation. With this understanding, we will follow the conventions of [Fj, and most of the 
properties of NBV([0,T]; R) that we use can be found there. 

It is easy to check that NBV([0, T];R) is closed under pointwise products and sums. 
Functions in NBV([0, T}; R) have bounded images, at most countably many discontinuities, 
and well-defined left-hand limits at each discontinuity. In particular, they are examples of 
cadlag functions {continue a droite, limites a gauche). (In fact, NBV([0, T]; R) is exactly the 
set of cddldg functions of bounded variation.) For any u G NBV([0, T]; R) and any s G [0, T], 
we set 

u(s-) = limit (£), 

v ' t/s W 

and we define the jump of u at s to be 

A s (it) = it(s) — u(s—). 

Note that u(0—) = and Aq(u) = u(0) by our conventions. 

It follows from elementary measure theory that for each u G NBV([0, T]; R), there is a 
unique signed Borel measure du on [0, T] characterized by 

du((o,6]) =u{b) -u(a), te[0,T\. 

Because this measure has atoms exactly at points t G [0, T] where u is discontinuous, we 
have to be careful to indicate whether endpoints are included or excluded in integrals. For 
example, we have the following versions of the fundamental theorem of calculus for a, b G 
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[0,T]: 



du = uCb) — u(a); / du — u(b) — u(a— ); 

(a,b] J[a,b] 



du = u(b—) — u(a); / du — u(b— ) — u(a— ). 

(a.,6) J[a,b) 

The tota/ variation of n, denoted by is given by either of two formulas: 

|| = sup |^ \u(xi) - u(xi- X )\ : = x < xi < ■ ■ ■ < x k = T j 

= / \du\. 

J[0,T] 

It follows from our conventions that ||m||oo < \\du\\. 

For u G NBV([0, T]; R), we will use the notation m_ to denote the function u_(£) = u(t—). 
Note that w_ has bounded variation, but is left-continuous rather than right-continuous. 

Lemma 2.3. For any u, v G NBV([0, T];R) and a, 6 G [0, T], £ae following integration by 
parts formula holds: 



udv+ V- du = u(b)v(b) — u(a)v(a). (2-11) 

(a,b] J \a,b] 

Proof. This follows as in (PJ Thm. 3.36] by applying Fubini's theorem to the integral du x 
dv, where Q is the triangle {(s,t) : a < s < t < b}. □ 

Lemma 2.4. Tne following product rules hold for m,d6 NBV([0, T]; R): 

d(nv) = udv + V- du 
= u~ dv + v du 

= udv + vdu-^2 A Sl (n) A Si (w)5 Si , 
j 

where 5 Si is the Dirac mass at Si, and the sum is over the countably many points Sj G [0, T] 
at which both u and v are discontinous. 

Proof. The first two formulas follow immediately from (12.111) and the definition of d{uv). 
For the third, we just note that the measure (v — vJ)du is supported on the set of points 
where u and v are both discontinuous, and for each such point Sj, 

(v(si) - v-(si))du({si}) = (v(si) - v(si-))(u(si) - w(s— )) 

= A Si (u)A Si (v). 

□ 

We will be interested primarily in vector- valued functions. We let NBV([0, T]; R n ) denote 
the set of functions v: [0, T] — > R" each of whose component functions is in NBV([0, T]; R), 
and NBVQO, T];M) C NBV([0, T]; R n ) the subset of functions taking their values in M. 
The considerations above apply equally well to such vector-valued functions, with obvious 
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trivial modifications in notation. For example, if v, w £ NBV([0, T]; R n ), we consider dv 
and dw as R n - valued measures, and Lemma 12.31 implies that 



If 7 G NBV([0, T];M), we say a function v G NBV([0, T] ; R n ) is a vector /ieW a/on# 7 if 
v(£) G Tj(t)M for each i G [0, T]. This is equivalent to the equation (v(t), n(7(£))) = for all 
t, or more succinctly (v, n o 7) =0. Note that the fact that 7 takes its values in a bounded 
set, on which n is uniformly Lipschitz, guarantees that n o 7 G NBV([0, T]; R n ). 

We generalize the notion of covariant derivative for NBV vector fields by defining 

Dv — dv + (v_, d(n o 7))n o 7. 

One motivation for this definition is provided by the following lemma, which says that if 
v(0) is tangent to M and Dv is tangent to M on all of [0, T], then v stays tangent to M. 

Lemma 2.5. Suppose 7 G NBV([0,T]; M) and v G NBV([0,T]; R n ). J/v(0) G T 7(0) M and 
(Dv, n o 7) = 0, then v(t) G 7^ (t) M /or all t G [0, T] . 

Proof. Using Lemma [2 A\ we compute 



2.3. An Existence and Uniqueness Theorem. The main purpose of this section is to 
prove the following theorem. 

Theorem 2.6. Let M C R n be a smooth, properly embedded hyper surf ace, and let 7 G 
NBV([0,T]; M). For any v G T^M, there exists a unique NBV vector field v along 7 that 
is a solution to the following (measure-valued) ODE initial-value problem: 



Before proving the theorem, we will establish some important preliminary results. We 
begin by dispensing with the uniqueness question. 

Lemma 2.7. Let 7 G NBV([0, T\ \ M). I/v,v G NBV([0, T]; R") are both solutions to fl2~T2l) 
with the same initial condition, they are equal. 

Proof. Suppose v is any solution to (12.121) . Observe that Lemma 12.51 implies that v(i) is 
tangent to M for all t, so (v, n o 7) =0. Let R = ||t lloo ? so that 7 takes its values in Mr. 




= (Dv, n o 7} 
= (dv, n o 7) + {^(v_, d(n o 7))n o 7, n o 7^ 
= d(v, n o 7) — (v_, d(n o 7)) + (v_, d(n o 7))(n o 7, n o 7) 
= d(v, n o 7). 

Thus if (v(0), 11(7(0))) = 0, we find by integration that {v(t), n(j(t))) = for all t. 



□ 



Dv = (So^)vdt, 
v(0) = v . 



(2.12) 
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With K chosen as in Lemma EH define / G NBV([0,T];R) by f(t) = e" 2i "|v(t)| 2 . Then 
Lemma 12.41 yields 

df = e~ 2Kt (-2K\v\ 2 dt + 2(v, dv) - £(A s ,v, A Si v)5 s ) 
= e~ 2Kt (^-2K\v\ 2 dt - 2(v, n o 7 }(v_, d(n o 7)) 
+ 2(v, (S o 7 )v>dt - £(A Si v, A Si v)<^ 

= e" 2 ^ ^(<v, (5 o 7 )v) - i^|v| 2 )^ - J] I A Sl v| 2 ^ . 

Since (12. 4p shows that (v, (S o 7 )v) < i^|v| 2 , this last expression is a nonpositive measure 
on [0, T]. Integrating, we conclude that f(t) < /(0), or 

|v(t)| 2 <e 2 ^|vo| 2 - 

In particular, the only solution with initial condition v = is the zero solution. Because 
(I2.12p is linear in v, this suffices. □ 

To prove existence, we will work first with finite approximations. Define a finite trajectory 
in M to be a function 7 G NBV([0, T];M) that takes on only finitely many values. This 
means that there exists a partition {0 = t < t\ < ■ ■ ■ < t m = T} of [0, T] such that 7 
is constant on for each z. For such a function, c?7 = J2iLo At 4 ( 7 )<5t t and ||g? 7 || = 

E£o 1^(7)1- 

Suppose 7 is a finite trajectory in M and v G 7^( )M. Let = to < • • ■ < t m = T be a 
finite partition of [0,T] including all of the discontinuities of 7, and write Xi = 7 (£j). Define 
v: [0,T] -> R n by 

v(t) = e^^^TT^e^-^-^^fc-i^^j ■ ■ .e fe ~^i7r Xl e (tl_to) *»°Vo, (2.13) 

where k is the largest index such that < £. Observe that the definition of v is unchanged 
if we insert more times in the partition. 

Lemma 2.8. Let 7 : [0, T] — > M be a finite trajectory. For and any vo G T^^M , the map v 
defined by (12.131) is the unique solution to (12.121) . and satisfies 

|v(t)|<e ct |v |, (2.14) 
||dv|| < C, (2.15) 
where C is a constant depending only on M , T , and 
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Proof. An easy computation shows that 

m 

dv = {S o 7 ) v dt + Yl - v(t<-)) ^ 

1=0 
m 

= (5 o 7 )v tft + £(v(ti-), n( 7 (U)) - n( T (t-)))n( 7 (t i ))5 ti , 

8=0 

from which it follows that v solves (I2.12p . 

To estimate |v(t)|, observe first that the operator norm of each projection ir x is equal to 
one. Let K be the constant of Lemma \2. II for R = \\dj\\. Using (I2.6p . we have the following 
operator norm estimate for any finite collection of points xi, . . . , Xj G Mr and real numbers 
h,...,lj G[0,T]: 

\\e l i Sx i o7i Xj o---o e hSx i o ti Xi || < e Kl > ■ ■ ■ e Kh = e K (h+-+h)^ ( 2 .16) 
Applying this to the definition of v proves (I2.14p . Then, using (I2.14p and fl2.10p . we estimate 



„ m 

iidvii= / \(soj)v\dt+Y: 



v{t~), n( 7 (f»)) - n( 7 (t,-)))n( 7 (^)) 



„ m 

< / Ke Kt dt + J2 e^Kl^U) - 7 (t-) 



'10,T] i=0 
< C(l + ||d 7 ||). 

□ 

Lemma 2.9. Suppose 7 and 7 are any finite trajectories in M defined on [0, T] and starting 
at the same point, and v ; v are tae corresponding solutions to (12.121) . There is a constant 
C depending only on M, T, || 7 ||oo; ond || 7 ||oo such that the following estimate holds: 

||v - v||oo < C (1 + ||d 7 || + ||d 7 ||) || 7 - 7||oo|Vo|. 

Proof. Lemma 12.81 shows that \\v\\oo an d ||v||oo are both bounded by C|v | for some C 
depending only on M, T, || 7 ||oo, and || 7 ||oo- Fix t e [0,T], and let = to < ■ • • < tk < t 
denote a finite partition that includes all of the discontinuities of 7 and 7 in [0,t]. We 
introduce the following shorthand notations: 

Xi = j{ti), Xi = 7 (tj), 

Si = S{Xi), Si = S(xi), 

Xi = Kxi: Ki = 7T~. 

Observing that 7r vo = v and 7^+1 v(t) = v(i), we can write v(i) — v(i) as a telescoping 
sum: 

v(t) - v(t) = e lkSk n k ■ ■ ■ e h+A+1 n i+1 (e liSi 7n - K i+1 e hSl \ 7f * - - - e illSl 7fie io<So v . 

i=0 ^ ' 
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By (I2.16p . the compositions of operators before and after the parentheses in the summation 
above are uniformly bounded in operator norm by e KT . Therefore, 

k 



|v(t)-v(t)|<e 2 ^E 



i=0 



vr-i+i o i e hSl o m - 7c i+1 o e kSl | O 77, 



l v o|- 



Using the fact that Si and 7Tj commute, as do Si and 7fj, we decompose the middle factors 
as follows: 

n i+1 o (e USi o 7Tj - Tr i+1 o e li ^ o rr { = ir i+1 o tt { o (e liSi - e li ^ o ^ 

+ %io( 7 r i -i !+1 )o^oe' ,Si . 

We will deal with each of these terms separately. 
For the first term, (I2.9I) implies 

and after summing over i, we find that this is bounded by KT^ — 7||oo. For the second 
term, Lemma 12.21 allows us to conclude that 

7T;+1 O (TTj - TTi+i) O 7fj O e' I<Sl 



< Ce XT ||7 - 7||oo (\x i+1 - Xi| + \x i+1 -Xi\). 



After summing, this is bounded by Ce KT ^ r y — 7 



+ 



. This completes the proof. 

□ 



Lemma 2.10. Let 7 6 NBV([0, T];M) 6e arbitrary. For any e > 0, there exists a finite 
trajectory 7: [0, T] — > M stzc/i that H7 — 7H00 < £ ond ||<i7|| < 



Proof. Let £ be given. Since 7 is cadlag, for each a 6 [0, T], there exists 5 > such that for 

f e [0,21, 

te[a,a + 5) =^ \j(t) -7(a) I < e, (2.17) 

t€(o-i,o) =*► l7(*)-7(OI<|- (2-18) 
By compactness, we can choose finitely many points = Oq < ai < . . . a m = T and 
corresponding positive numbers So,...,S m so that [0, T] is covered by the intervals (a^ — 
5j, aj + 5i), i = 1, . . . , m. Because they are a cover, for each i — 1, . . . , m we can choose 6j 
such that 

6, G (oj-i, + Si-i) n (a* - Si, ai). 
Now define a finite trajectory 7 : [0, T] — > M by 

7(ai_i), t e [oi_i,6i), 



7(0 



7(6i), t E [bi,ai). 
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It is clear from the definition of the total variation that ||d7|| < ||d7||. We will show that 
||7 - 7l|oo < e. 

Let t G [0, T] be arbitrary. For some i, either t G [aj__i,6j) or t G [6^, aj). In the first case, 
since [aj_i,6j) C [aj_i,aj_i + <5j_i) by construction, (I2.17P yields 

|7(t)-7(*)l = l7(*)-7(ai-i)l<e- 
On the other hand, if t G cij) C (a, — 5j, dj), (12.181) yields 

|7(*)-7(*)l = l7(*) -7(MI 

< | T (t)- T (a-)| + l7K-)-7(6i)l 

< 2 + 2' 

so we reach the same conclusion. □ 

Lemma 2.11. For any 7 G NBV([0, T]; M), there exists a sequence of finite trajectories 
7W : [0, T] — > M satisfying Wdj^W < \\d"y\\ and converging uniformly to 7. 

Proof. This is an immediate consequence of Lemma 12.101 □ 

Now we can prove the existence and uniqueness theorem. 

Proof of Theorem \2.b\ Given 7 as in the statement of the theorem, let 7^ be a sequence of 
finite trajectories converging uniformly to 7 as guaranteed by Lemma 12.111 For each k, let 
v( fc ) be the solution to (12.121) for 7 = 7W, as defined by (I2.13p . Then Lemma T2.9I guarantees 
that the sequence v( fc ) is uniformly Cauchy, and hence there is a limit function v : [0, T] — > R n 
such that v( fc ) — > v uniformly. It is straightforward to check that v G NBV([0, T]; R n ). 
Moreover, since each is tangent to M and v( fe ) — * v uniformly, it follows that v is also 
tangent to M. 

We need to show that v solves (12 . 1 2[) for 7. It suffices to show for any w G NBV([0, T]; R n ) 
that 

/ (w, dv) = — (w, n o 7)(v_, d(n o 7)) + / (w,(S oj)v) dt. 

J[0,T] J[0,T] ' J[0,T] 

If we write w = w T + w^, where w T is tangent to M and w 1 - is orthogonal to M, this is 
equivalent to the following two equations: 

[ ( w ± ,rfv) = -/ (w ± ,no 7 )(v_,d(no 7 )), (2.19) 

J[0,T] J[0,T] 

f (w T ,rfv)= f (w T ,(5o 7 )v)cit. (2.20) 

J[0,T] J[0,T] 

Because w -1 is proportional to n, w x = (w -1 , n o 7)11 o 7. The fact that v is tangent to M 
means that (n o 7, v) = 0, from which we conclude 

= d(n o 7, v) = (n o 7, dv) + (v_, d(n o 7)). 
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Therefore, 

(w^, dv) = ^(w^, n o 7)11 o 7, dv^ 
= (w -1 , n o 7)(n o 7, dv) 
= -(w ± ,no 7 )(v_,rf(no7)), 

from which (12 . 19[) follows. 

On the other hand, from Lemma 12.31 we conclude that 

/ (w T ,rfv) = (w T (T),v(T)> - / (v_,dw T ) 

J[0,T] J[0,T] 

= lim [ (w T (T),v( fc )(T)) - / (vi fe) ,rfw T ) ) 

k^oo y J[0,T] J 

= lim / (w T ,dv (fc) ) 

A; -+00 J[0,T] 

= lim f- / (w T ,no 7 W)(v (fc) ,d(no 7 ( fc ))) 

fe-too^ 7[0,T] 

+ / (w T ,(So 7 ( fc ))v«)cA 

y[o,T] y 

Since (w T , 1107^)} converges uniformly to (w T , 1107) = 0, and the measures (v_ , d(no'-f^)) 
have uniformly bounded total variation, the first term above vanishes in the limit. Since both 
v( fc ) and 5 o 7'^ converge uniformly, the last term above converges to /[ 0T ](w T , (S °7)v) dt. 
This proves fl2T20|) . □ 

2.4. Stability. In this section, we wish to address the stability of the solution to (12.121) 
under perturbations of the trajectory 7. For applications to probability, we will need to 
consider perturbations in a weaker topology than the uniform one. 

We define a metric d$ on NBV([0,T]; R n ), called the Skhorokhod metric, by 

^5(7,7) = pf max(||7 -70 A||oo, ||A - Id H^) , 

where A is the set of increasing homeomorphisms A: [0, T] — > [0,T]. We wish to show that 
the solution to (I2.12p is continuous in the Skorokhod metric, as long as we stay within a set 
of trajectories with uniformly bounded total variation. 

Because the Skorokhod metric is not homogeneous with respect to constant multiples, it 
will not be possible to bound d$(v, v) directly in terms of 0^(7,7). For this reason, we will 
work instead with the solution operator, for any NBV trajectory 7: [0,T] — > M, this is the 
endomorphism- valued function A: [0, T] — > End(R n ) defined by 

A(t)v = v(t), 

where v is the solution to (I2.12p with initial value v , and extended to an endomorphism 
of R n by declaring ^4(t)n 7 ( ) = 0. As before, ||-4(t)|| will denote the operator norm of A(t), 
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and we set 

M|| oo = sup{M(t)||:tG[0,T]} 

= sup jl^jl : t G [0,T], v G T y{0) M, v ^ o| . 

It follows easily from the results of the preceding section that for any 7 G NBV([0,T]; M), 
the solution operator A is in NBV([0, T]; End(R n )), and Lemma T2.9I translates immediately 
into the following estimate. 

Lemma 2.12. Suppose 7 and 7 are any finite trajectories in M defined on [0, T] and starting 
at the same point, and A, A are the corresponding solution operators. There is a constant 
C depending only on M , T, ||g?7|| ; and \\drf\\ such that the following estimate holds: 

\\A-A\loo < CH7-7IU 

Next we need to examine the effect of a reparametrization on the solution associated with 
a finite trajectory. 

Lemma 2.13. Lei 7: [0, T] — > M be a finite trajectory, let A: [0,T] — > [0, T] be an increasing 
homeomorphism, and let 7 = 7 o A. There is a constant C depending only on M, T, and 
WdjW such that the solutions v and v to (12.121) associated to 7 and 7 with the same initial 
value v satisfy 

||v-vo Alloc < C||A - Id ||oo|vo|- (2.21) 



Proof. As in the proof of Lemma [2.91 fix t G [0, T] and let = to < t\ < ■ ■ ■ < tk < t be the 
points in [0, £] at which 7 is discontinuous. Set t^+i = t 1 x i = 7^), and = A(£,), so that 7 
and 7 are given by 

7(t) = Xi if ^ < t < t;+i, 

7(t) = ^ if tj < t < t ;+ i. 

We will also use the notations 
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v(A(t)) - v(t) = (ld-e ( * fe + 1 - ?fc + l)Sfc ) v(A(t)) 
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+ E 



1=1 

e 



e K K 7T fc • • ■ e 1+1 l+1 7T i+1 



TTi-ie 



ij_2<Si_2 



7r, ; e 



e' lSl 7rie ioiSo v . 



By virtue of (12.81) . the first term is bounded by a constant multiple of |ifc+i — tfc+i| |v | < 
|| A— Id ||oo| v o|- As before, the compositions before and after the parentheses in the summation 
are uniformly bounded in operator norm, so we need only estimate the sum 
k 



E 

i=l 



e (t i+1 ~ti)Si Q Q Ji-xSi-x _ e hSi Q Q g^-i.-i)^-! 



Using the fact that 7Tj commutes with iSj, we can rewrite the i-th term in this sum as 



e liSi o 7T,- o ( e (**-*;)Si 



g(*i — *i)5t— 1 j g'j-l<Si-l 



< 



From (12. 6p and (12. 9p , this last expression is bounded by C 
i, we conclude that this is bounded by C||A — Id || c 



ti tj 



\xi — Xi-i\. Summing over 

□ 



Lemma 2.14. Suppose 7,7: [0, T] — > M are finite trajectories starting at the same point, 
and let A, A be the corresponding solution operators. There exists a constant C depending 
only on M, T, ||<i7||, and ||g?7|| such that 

d s (A,A) < 6^(7,7). (2-22) 

Proof. Let 5 = ds(j, 7) and let e > be arbitrary. By definition of the Skorokhod metric, 
there is an increasing homeomorphism A: [0,T] — > [0, T] such that H7 — 7 o A||oo < 5 + e 
and || A — Id || oo < 5 + e. Let A\ be the solution operator associated with 7 o A. Then 
\\A - AIL < C(5 + e) by Lemma [2ZCU and - A x o < C(<5 + e) by Lemma EE3J 
Thus by the triangle inequality, 

d s (A,A) < d s (A,Ai)+d s (Ai,A) 

< \\A- Ai\\oo + max (||.4 — Ai o A||oo, ||A - Id H^) 

< C(8 + e)+ max(C(5 + e), e). 

Letting e — > 0, we obtain 

d s (A,A)<2Cd s (rt,7). 

□ 



Here is our main stability result. 
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Theorem 2.15. Given positive constants R and T , there exists a constant C depending only 
on M, R, and T such that for any trajectories 7,7 G NBV([0, T];M) starting at the same 
point and with total variation bounded by R, the corresponding solution operators A and A 
satisfy 

d s (A,A) <Cd s (in). (2.23) 

Proof. By the argument in the proof of Theorem I2.6[ there exist sequences of finite trajec- 
tories converging uniformly to 7 and 7 whose solution operators converge uniformly to A 
and A, respectively. Thus for any e > 0, we can choose finite trajectories 7' and 7', with 
corresponding solution operators A' and A', such that 

II 7 tIIoo^^) II7 tIIoo^^) 

WA'-AWoo < e, H-4' -.A||oa < e. 

Then by the triangle inequality, 

d s (l\ l') < d s (j', 7) + d s (j, 7) + d g (7, 7') < d s (7, 7) + 2e. 
By Lemma 12 .14[ we have 

d s {A',A') < Cd s {irf) < Cd s (j,j) + 2Ce. 
Thus by the triangle inequality once more, 

d s (A, A) < d s (A, A') + d s (A', A') + d s (A', A) 
<e + {Cd s { 1 ^) + 2Ce)+e. 
Letting e — > completes the proof. □ 

2.5. Base trajectories of infinite variation. In the probabilistic context, we will have 
to analyze the situation when the base trajectory 7 does not have finite variation on finite 
intervals. We will now present an example showing that some of the results proved in 
this section do not extend to (all) functions 7 of infinite variation. Hence, arguments using 
piecewise-constant approximations in the probabilistic context will require some modification 
of our techniques. 

Example 2.16. Let M C R 2 be the parabola M = {(xx,x 2 ) G R 2 : x 2 = x 2 }, with the 

orientation of M chosen so that ||«S X || < 1 for all x G M. Let j(t) = (0, 0) for t G [0, 1], and 

for even integers j > 2, let 

f x, := for t G [2kj- 3 , (2k + l)r 3 ), k = 0, 1, . . . , f/2 - 1, 

lj(t) = Vi ■= H~\j~ 2 ), for t G [(2k + l)r 3 , (2k + 2)j- 3 ), k = 0, 1, . . . , j 3 /2 - 1, 
[ (r\r 2 ), fort = l. 

Clearly, jj — >■ 7 in the supremum norm on [0, 1], so ^(7^,7) — > 0. Let v = (1,0) and let 

Vj(t) be defined as in (12.131) . relative to jj. Similarly, let v(t) be defined by (I2.13P relative 

to 7. We have v(l) = e^ ' ) v ^ (0, 0). 
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There exists C\ > such that for all j > 2, z G T Xj M, we have |7T % -z| < (1 — C]j~ 2 )|z|, 
and similarly, \tt x .z\ < (1 — cij~ 2 )|z|, for z G T Vj M . This implies that for some C2 < 1, 
| v j (1)| = \{^xj o 7r yj ) :?3//2 vo| < c^. Hence, lim^oo Vj(l) = (0,0) ^ v(l). This shows that 
results such as Lemma 12.91 do not hold for (some) functions 7 which do not have bounded 
variation. 

3. Multiplicative functional for reflected Brownian motion 

Suppose D C R n , n > 2, is an open connected bounded set with C 2 boundary. Recall that 
n(x) denotes the unit inward normal vector at x G dD. Let B be standard rf-dimensional 
Brownian motion, G D, and consider the following Skorokhod equation, 

X t = x* + B t + C n{X s )dL s , for t > 0. (3.1) 
Jo 

Here L is the local time of X on <9.D. In other words, L is a non- decreasing continuous 
process which does not increase when X is in D, i.e., Jq 00 l£>(Xt)dLt = 0, a.s. Equation (13.11) 
has a unique pathwise solution (X, L) such that X t & D for alH > (see |LS] ). 

We need an extra "cemetery point" A outside R n , so that we can send processes killed at a 
finite time to A. Excursions of X from dD will be denoted e or e s , i.e., if s < u, X s , X u G dD, 
and X t £ dD for t G (s, u) then e s = {e s (t) = X t+S , t G [0, w — s)}. Let C( e s) = u — s be the 
lifetime of e s . By convention, e s (t) = A for t > (, so e t = A if inf{s > t : X s G dD} = t. 

Let a be the inverse local time, i.e., at = inf{s > : L s > t}, and S r = {e s : s < 
cr r }. Fix some r, e > and let {e 4l , e f2 , . . . , e fm } be the set of all excursions e <E £ r with 
|e(0) — e(£— )| > e. We assume that excursions are labeled so that tk < tk+i for all k and we 
let £ k — L tk for k — 1, . . . , m. We also let t = inf{t > : X t G <9-D}, ^0 = 0, l m +i = r, and 
A4 = 4 + i - 4. Let x k = e tk ((-) for A; = 1, . . . , m, and let x = X to . 

In this section, the boundary of D will play the role of the hypersurface M, i.e., M = dD. 
Recall that S denotes the shape operator and tt x is the orthogonal projection on the tangent 
space T x dD, for x G dD. For v G R n , let 

v rj£ = exp(A£ m S(x m ))n Xm ■ ■ ■ exp(A45(si))7r :ri exp(A£ «S(x ))7r a;o Vo. (3.2) 

Let A r:e be a linear mapping defined by v r £ = ^4 r e v . 

We point out that the "multiplicative functional" At discussed in the the Introduction is 
not the same as A r defined in this section. Intuitively speaking, A r = A ar , although we 
have not defined At in a formal way. 

Suppose that dD contains n non-degenerate (n — l)-dimensional spheres, such that vectors 
perpendicular to these spheres are orthogonal to each other. If the trajectory {Xt, < t < r} 
visits the n spheres and no other part of dD, then it is easy to see that A r)£ = for small 
e > 0. To avoid this uninteresting situation, we impose the following assumption on D. 
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Assumption 3.1. For every x G dD, the (n — 1)- dimensional surface area measure of 
{y G dD : (n(y),n(x)) = 0} is zero. 

Theorem 3.2. Suppose that Assumption \3.1\ holds. With probability 1, for every r > 0, the 
limit Ar '■= lim e ^o^r,e exists and it is a linear mapping of rank n — 1. For any Vo, with 
probability 1, ^4 r>£ v — > A r v uniformly on compact sets. 

Remark 3.3. Intuitively speaking, A r vo represents the solution to the following ODE, similar 
to (I2.12p . Let j(t) = X(a t ), and suppose that v G R n . Consider the following ODE, 



Then A T is defined by v(r) = A t vq. We cannot use Theorem 12.61 to justify this definition of 
A r because 7 ^ NBV([0, r]; dD). See [X], [IK1] or [H] for various versions of the above claim 
with rigorous proofs. Those papers also contain proofs of the fact that A r is a multiplicative 
functional of reflected Brownian motion. This last claim follows directly from our definition 

of Ar. 

Remark 3.4. Recall that B is standard <i-dimensional Brownian motion and consider the 
following stochastic flow, 



where L x is the local time of X x on dD. The results in |LSj are deterministic in na- 
ture, so with probability 1, for all x G D simultaneously, (13.31) has a unique pathwise 
solution (X X ,L X ). In a forthcoming paper, we will prove that for every r > 0, a.s., 



The rest of this section is devoted to the proof of Theorem 13.21 We precede the actual 
proof with a short review of the excursion theory. See, e.g., [M] for the foundations of the 
theory in the abstract setting and [BuJ for the special case of excursions of Brownian motion. 
Although [Buj does not discuss reflected Brownian motion, all results we need from that book 
readily apply in the present context. 

An "exit system" for excursions of the reflected Brownian motion X from dD is a pair 
(L£, H x ) consisting of a positive continuous additive functional L* t and a family of "excursion 
laws" {H x } x( zdD- In fact, L* t = L t ; see, e.g., |BCJj . Recall that A denotes the "cemetery" 
point outside R n and let C be the space of all functions / : [0, 00) — > R n U {A} which are 
continuous and take values in R n on some interval [0, £), and are equal to A on [£,00). 
For x G 3D, the excursion law H x is a a-finite (positive) measure on C, such that the 
canonical process is strong Markov on (t , 00), for every t > 0, with transition probabilities 
of Brownian motion killed upon hitting dD. Moreover, H x gives zero mass to paths which do 
not start from x. We will be concerned only with "standard" excursion laws; see Definition 



Dv = (Soj)vdt, 



V(0) = 7T X . V . 




for t > 0, 



(3.3) 



lirn^o sup^^ (X^ +£V - X%)/e - Av = 0. 
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3.2 of [Buj . For every x G dD there exists a unique standard excursion law H x in D, up to 
a multiplicative constant. 

Recall that excursions of X from dD are denoted e or e s , i.e., if s < u, X S ,X U G dD, 
and X t £ dD for t G then e s = {e s (t) = X t+S , t G [0,w — s)} and C( e s) = u — s. By 

convention, e s (t) = A for t > £, so e t = A if infjs > t : X s G <9.D} = 1 

Recall that at = inf{s > : L s > t} and let / be the set of left endpoints of all connected 
components of (0, oo) \ {t > : X t G dD}. The following is a special case of the exit system 
formula of [M], 



E 



E Vt ■ f(e t ) 



E 



/ V a H x ^\f)ds = E V t H Xt (f)dL t , (3.4) 
Jo Jo 



where V t is a predictable process and / : C — > [0, cxd) is a universally measurable function 
which vanishes on excursions et identically equal to A. Here and elsewhere H x (f) = J c fdH x . 

The normalization of the exit system is somewhat arbitrary, for example, if (L t ,H x ) is 
an exit system and c G (0, oo) is a constant then (cL t , (l/c)H x ) is also an exit system. Let 
Pi) denote the distribution of Brownian motion starting from y and killed upon exiting D. 
Theorem 7.2 of [Buj shows how to choose a "canonical" exit system; that theorem is stated 
for the usual planar Brownian motion but it is easy to check that both the statement and 
the proof apply to the reflected Brownian motion in R n . According to that result, we can 
take LI to be the continuous additive functional whose Revuz measure is a constant multiple 
of the surface area measure on dD and H x, s to be standard excursion laws normalized so 
that 

H x (A)=\^Pl +S ^\A), (3.5) 

for any event A in a a-field generated by the process on an interval [t , oo), for any t > 0. 
The Revuz measure of L is the measure dx/(2\D\) on dD, i.e., if the initial distribution 
of X is the uniform probability measure \i in D then E^^ 1 l A (X s )dL s = f A dx/(2\D\) for 
any Borel set A C dD, see Example 5.2.2 of [FOTj. It has been shown in |BCJj that 
(Lj , H x ) = (L t , H x ) is an exit system for X in D, assuming the above normalization. 

Proof of Theorem \3.2 . The overall structure of our argument will be similar to that in the 
proof of Lemma 12.91 

We will first consider the case r — 1. Let Ej = 2 _J , for j > 1. Fix some j for now and 
suppose that e' G [ej + i,ej). Let 

{ e t{> e 4> • • • > } = {e G : |e(0) - e(C~)| > 
{e t[ ,ee 2 , e^} = {e G £ x : |e(0) - e(C-)| > e'}. 



m 



j ■ 



We label the excursions so that t J k < t J k+l for all k and we let l\ = L tJ for k — 1, 
Similarly, ^ < t^+i f° r all ^ an d i' k — L t i for = 1, . . . , m'. We also let t J = t' = inf{t > 
: X t G dD}, ^ = ^o = 0, = C'+i = X > A£ i = 4+1 " 4. and A 4 = 4+i ~ 4- Let 
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e t j ((—) for k = 1, . . . ,rrij, and x' k = e t > ((—) for k = 1, . . . ,m'. Let Xq = X t j, and 



Let 7 J (s) = 4 for s G [£{, £{ +1 ) and fc = 0, 1, ... , m h and 7 J (1) = 7 j (4J- Let f (s) 



,1 ; 



for s G [4, f fc+1 ) and fc = 0, 1, ... , m' , and = 7 ' (CO- 
For v G R", let 

= exp(A4 5(< ^ • ■ ■ exp(Afi5(xi))7r , exp(A^(^)K f) v , 
v' = exp(A£' m ,S(x' m ,))ir x > ■ ■ ■exp(A£[S(x' 1 ))Ti x ' exp(A£^S(x' ))'K x ' v . 

Let = £ < " " " < Art+i = 1 denote the ordered set of all f k s, < k < rrij + 1, and ^.'s, 
< /c < m' + 1. In the definition of £^s, we followed the proof of Lemma [2.91 word by word, 
for conceptual consistency, although the set of £k's is the same as the set of £' k s. 

We introduce the following shorthand notations, A» = £i + \ — £i, 

Xi = r y\£i), Xi = j'(£i), 

Si = S[Xi), Si = S[Xi), 

Observing that 7To7fo v o = 7fo v o and 7r m +iv' = v', we can write v- 7 — v' as a telescoping 
sum: 



\ j - v' = e Am5m 7r m • • •e iii+l0i+1 7T i+ i [ e" ,c, *7ri - t i+1 e"' lV% )*,■■■< l » •f-^"T v . 

i=0 



3 Ai5i~ o A 5 ; 



By (I2.16p . the compositions of operators before and after the parentheses in the summa- 
tion above are uniformly bounded in operator norm by a constant. Therefore, for some c\ 
depending only on D, 



m 

|v J - v'| < Ci^ 

i=0 



7Tj+l ° I e AA OTli- 7l i+1 O e A ' S * I O 77, 



l v o|- 



(3.6) 



Using the fact that Si and 7Tj commute, as do Si and TTj, we decompose the middle factors 
as follows: 



7Ti+l O I e AA O TTi - 7T i+1 O , S * ) O 7T; = 7T i+1 O 7T; O f ( A ^ - < A - S | O 7T ; 



(3.7) 



+ 7T i+ i O (7Tj - 7T i+ l) O 7Ti O e 

We will deal with each of these terms separately. 
For the first term, we have by (12.91) . 

7r i+1 o m o (e AlSl - e A '^ o 7fj < e Al,Sl -e Al ^ <c 2 A i |a; i -£ i 

For the second term, Lemma 12.21 and (I2.6P allow us to conclude that 

VTj+l O (m - TTi+i) O TTj O e AllSl 



(3.f 



3^ | 






+ 




— 




| 






+ 




— x i+l 





,Aj5< 
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We will now estimate Esup 0<i<m \xi — Xi\. Suppose that Xi ^ Xi for some i. Then there 
exist k\ and k 2 such that £ k < £' k2 < £ J kl+ i, Xi = x ki , and Xi = x' k . Hence, 

{\xi — Xi\ > a} C [J < sup \x k — X t \ > a > . (3.10) 

k l*i.+CK)<*<*UiA e9D J 
Intuitively speaking, the last condition means that the process X deviates by more than a 

units from x k (the right endpoint of an excursion ej ), when X is on the boundary of D, 

at some time between the lifetime of this excursion and the start of the next excursion in 

this family, e.j 

Since dD is C 2 , standards estimates (see, e.g., |Buj ) show that for some a ,c 5 > 0, all 
x G dD and a G (0, ao), 

l/(c 5 a) < H x (|e(C— ) -x\> a) < c 5 /a. (3.11) 

It follows from this and (13.41) that there exists cq so large that for any stopping time T and 
a G (0,a ), 

P (3e s : |e s (C-) - e s (0)| > a, s G (T, a[L T + c 6 a))) > 3/4. (3.12) 
Let Tjs< x ,a) be the exit time of X from the ball B(x, a) in R n with center x and radius a. 
Routine estimates show that for some C7, a x > 0, and all a G (0, a x ) and x G dD, 

P x {L{r B(X}Cra) ) > c 6 a) > 3/4. (3.13) 

Let T k,o = tl and 

T{ i+1 = mf{t > Tl- : X(t) G dD, \X(t) - X{TQ\ > c 7 e 3 }, 1 > 0. 

According to (I3.13p . the amount of local time generated on (T k ,T kl ) will be greater than 
c^Ej with probability greater than 3/4. This and (13.121) imply that there exists an excursion 
e s with \e s ((— ) — e s (0) | > Ej and s G (T k ,T k l ), with probability greater than 1/2. By the 
strong Markov property, if there does not exist an excursion e s with \e s ((—) — e s (0) | > £j 
and s G (T k0 ,T ki ) then there exists an excursion e s with |e s (£— ) — e s (0)| > Ej and s G 
{T ki ,T ki+1 ), with probability greater than 1/2. Let M k be the smallest % with the property 
that there exists an excursion e s with \e s ((—) — e s (0) | > Sj and s G {T ki ,T ki+1 ). We see 
that Ml is majorized by a geometric random variable Ml with mean 2. Note that 

\X(Tl i+1 ) - X{TQ\ < (c 7 + 1)6, = ess,, 

for i < Ml. Therefore, 

sup \x{- X t \ < csMlej. (3.14) 

tl+aei)<t<t{ +1 ,X t edD 

It is easy to see, using the strong Markov property at the stopping times t k , that we can 
assume that all {M 3 k , k > 0} are independent. 
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Consider an arbitrary (3i < —1 and let rij = Ej 1 . For some c 9 > 0, not depending on j, 
P (max c,M{e 3 > W) = 1 - (1 - (1/2)*)"' < 1 1 . ,^ * \^ ^\ (3-15) 



CgUj (1/2)* iii>(3 1 j. 



Let p be the diameter of D and ji be the largest integer smaller than logp - By (I3.15p . 
for any /3 2 < 1, some Ci 2 < oo, and all j > jx, 



E ( max cgMlej J < E max cgMlej 



< c & i£ j+ J2 c 8 ie j c 9 n j (l/2y 

i</3ij i>/3ij 

< c 10 Sj(\ogej) 2 + cnSjl loge^ < ci 2 ef . (3.16) 

Let N e be the number of excursions e s with s < o\ and |e s (0) — e s ((—)\ > e. For e = Ej, 
N e = rrij. Then (13 .4[) and (13.111) imply that N £ is stochastically majorized by a Poisson 
random variable N e with mean Cis/e, where c\3 < oo does not depend on e > 0. We have 
Eexp(iV e ) = exp(ci 3 e _1 (e — 1)), so for any a > 0, 

P{N £ >a)< P{N £ >a) = P(exp(iV £ ) > exp(o)) < exp^e -1 - a). 

Standard calculations yield the following estimates. For any /3 3 < —1, /? 4 < 0, 8\ > 0, some 
<5 2 G (0, and all 83,84 £ (0,5 2 ), 

P(N S3 > 8§*) < 81 (3.17) 

and 

CSl E ( m {^})- sl (318) 

It follows from (I3.14p . (13.161) and (13.171) that, for any /3 2 < 1, some Ci 6 , and j > ji, 



E ( max sup \x 3 k — X 4 | 

"^<™J t i+«e{)<t<t{ +1 ,X t edD 



< E I max sup \x k — X t \ + p P( m j > ra j) 

V - fe - n, '*i+c(eixt<4 +1 ,x t6 sD / 

^ E (o^- C8M ^) +Cl5£ ? 

<c 12 ef +c 15 4<c 16 ef. (3.19) 
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Note that ET=o A i = 1 - This > (EH, (M and f[3T0D imply that, 
E( sup £ 7r i+1 o7r i0 fe A ^ ! - e A ^) o ^ j < E ( sup ]T caA^a* - Z;| j (3.20) 



\ej+i<e'<£j i=0 



< E sup max \xi — Xj| V"c 2 Aj = c 2 E sup max |xj — x/\ 

Ve J+1 < e '< ej °< i < m i=0 / V^+i<^'<^ °- i - m / 

< c 2 E max sup \xj, — X t \\ < ciqSj. 

\°^ m i i+«4)<t<ti +1 ,x t e8D " J 

We will now estimate the right hand side of (I3.9p . We start with an observation similar 
to (I3.10p . Suppose that Xi ^ Xi + \ for some i. Then there exists k\ such that Xi = x 3 ^, and 
x i+i = 4i+i- Note that fci's corresponding to distinct i's are distinct. Hence, 

{\xi - x i+l \ > a} (3.21) 



C |J \ SU P 14 - x t\ > a/2 > U (|e tJ (0) - ((-){ > a/2) 

k [t{+<;(el)<t<t{ +1 ,x t edD J L k+1 k+1 > 

C (J J sup \4 - X t \ > a/2 \ U |J {|e^+i(0) - e fJ+ i(C-)| > a/2} . 

k {t{+t;(e{)<t<t{ +1 ,x t edD J fe k k 

Similarly, suppose that Xi 7^ x i+ i for some i. Then there exists fc 2 such that X{ = x' k2 , and 

= 4 2 +i- Again, fc 2 's corresponding to distinct z's are distinct. Hence, 

>a}c[J sup 14-^*1 > a / 2 MK +1 ( W +1 (C-)I > G / 2 }' 

Since £j + i < e' < Ej, this implies that, 

{\xi-x i+1 \>a} (3.22) 
C U SU P \4-X t \>a/2\u |J {|e ti+ i(0)-e ti+1 (C-)| > a/2\ . 

0<k< mj [tl+C(e J k )<t<tl +1 ,X t edD J 0<fc<m i+ i & 

It follows from (13~T0|) . (ET2T1) and (13T221 that 

sup {\ x i+l- x i\\ x i-Xi\ + \ x i+l-Xi+l\\%i+l-Xi\) {3.23) 

e J+ i<£'<£j 0<i< m 



2 



< 4 E L™ax. sup |4 -X t \ 

0<i<m y^ m J t{+C(e{)<t<t{ +1 ,X t &dD 

+ 8 max sup \x j k — X t \\ ^ le^+i (0) — e t j+i((— ) 

\0<k< mj f k+c{e J k)<t<t J k+i! x tedD J \o< k <m j+1 " k 



4(m + 1) max sup |4 — X-, 

\°- k - mj i+C(4)<t<ti +v x t edD 



2 



+ 8 oS sup |x * ~ Xt| £ l e ti+ l(0) ~~ e *i +l(C " 
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We have the following estimate, similar to (I3.16p . For any /3 5 < 2, some Ci 9 < oo, and 
3 > Ji> 



< J] (c 8 ie j ) 2 + J] (c 6 ie i ) 2 c & n J -(l/2) < 

i<Pij i>Pij 

< c l7 e)\ log ^| 3 + c 18 £|(log £i ) 2 < ci 9 ef . (3.24) 
We now proceed as in (I3.19p . For any /3 5 < 2 and j > j\, 

2 



E [ max sup \x 3 k — X t 

°- k - mj ti+C(ei)<t<ti +1 ,x t edD 



2 

■2i 



< E max sup \x 3 k — X t \ + p P(nij > rij) 

\°- k - ni i+C(ei)<t<ti +1 ,x t edD J 

<c 19 ef + c w e)<c 21 ef. (3.25) 

Recall that m is random and note that m < mj +1 . We obtain the following from (13.171) 
and (13.25j) . for any ft? < 1, by choosing appropriate (3$ < 2 and /3 6 < —1, 

2N 




E [ (m + 1) ( n max sup |4 - X t \\ (3.26) 

" ^K(ei)<*<**+i>**eM> / / 

\4-Xt\\ | +pgP(m+l>ef) 

" 1 - mi iP k +^ k )<t<A +1 ,x t edD J ) 

Next we estimate the second term on the right hand side of (13.231) as follows. The number 
of excursions e.j+i with |e,j+i(0) — e.j+i((— )| G [£j + i,£j] is bounded by m i+ i, so 

k k k 

j+i 

V |e^+i(0) - e^+i (C— ) I < m i £ i-i- 

k k 



Multiplicative functional for reflected Brownian motion 24 

Hence, for any j3 g < 0, we can choose /3 8 < 0, (3\ < — 1 and c 2 3 < oo so that for all j > j^, 
max sup |x£-X t | £ |ej+i (0) - ej+i(C-) 

V " " J ' ? k +C(ei)<t<i? k+1 ,X t e8D J Vo<fe<m J+1 

/ \ £i 

< max sup — X t | mjS^i 

Y- k - mi ti+<(4)<t<i +1 ^edD ' ) i= h 

I . \ 3+1 3+1 

\°- k - m > ti+ael)<t<ti +1 ,x t edD J i= h i=n l - 

< [max sup \4 - X t | j 2(j - j x )ef ^ + 2 Po £ m,l s e, 



- C23£ l 9 ( SU P 14 - ^1 ] + 2 Po I] m^l. s 8 ,£i. 



This, (13.191) and (I3.18P imply that for any j3 w < 1, by choosing an appropriate /3 2 < 1 and 
fa, < 0, we obtain for some c 2 6 < oo and j > ji, 

El max sup \4 - X t \) I ]T \e , +1 (0) - e. +1 (C-)| I (3.27) 

V - - * t> k +(H)<t<ti +1 ,x t edD J \o<k< mj+1 j 



< C2 3 ef E max sup |zj - X t | ] + 2p E I ^ ^i 1 /™^.^^; 

3+1 

< C24£j 9 £j 2 + C25 ^ e^-Si < C26£j 10 . 

»=j'l 

We combine (13.231) . (13.261) and (13.271) to see that for any /?io < 1, some C27 < 00 and all 

3 >ji, 

E SUp _ X i\ \ X i ~ %i\ + ki+l ~ ^+1 1 1^+1 _ ^1) - C 27 ef°. 

\£j+i<s'<£j 0<i<m J 

We use this estimate and (13.91) to see that for any j3 w < 1, some C27 < 00 and all j > ji, 

(3.28) 



'3+1 



E sup Y 

\£j + l<£'<£j i=Q 



n+\ o (-Ki - ir i+1 ) o m o e AlS ' 



<E SUp J^Ci (\x i+ i - Xi\ \Xi - Xi\ + \x i+ i - X i+1 \ \x i+1 - Xi\) < C 27 £ 

\ej+i<£'<£j i=Q / 



P10 
3 ' 



It follows (13. 6p . (13. 7p . (13.201) . and (I3.28P that for any (3 W < 1, some c 28 < 00 and all j > ji, 
E f sup \v j - V\) < c 28 ^ 10 |v | = c 28 2-' 3l «|v |. 
This implies that J2j>j 1 E (sup e _. +i<e / <£j . |v J — v'\j < 00, and, therefore, a.s., 

( sup |v J ' - v'| ) < 00. (3.29) 

3>3i \ e 3+i^ e ' <£ J / 
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We extend the notation v' from e' in the range [ej + i,Ej) to all e' > 0, in the obvious way. 
It is elementary to see that (13.291) implies that vi := lim e ^o v' exists. For every e' > 0, the 
mapping v — > v' is linear, so the same can be said about the mapping v — > vi := Alv . 

Note that the right hand side of (13. 6p corresponding to r G [0, 1) is less than or equal to 
the right hand side of (I3.6P in the case r — 1. Hence, we can strengthen (13.291) to the claim 
that a.s., 

sup sup | — v^.| < oo, 

where v J r and v' r are defined in a way analogous to v- 7 and v', relative to r G [0, 1]. The 
analogous argument shows that for any integer r > 0, a.s., 

sup sup \v J r — v' r \ I < oo. 

j>ji \ < r < r 0£j+i<£'<£j / 

We use the same argument as above to conclude that for any v , with probability 1, Ar >£ v — > 
A r Vo uniformly on compact sets. 

It remains to show that A r has rank n — 1. Without loss of generality, we will consider 
only r — 1. Recall definition (13.21) of v r and note that 7r xo v G T XQ dD. It will suffice to show 
that for any w G T XQ D such that w ^ 0, we have ^4iW ^ 0. 

Recall the definition of x-(.'s and related notation from the beginning of the proof. Recall 
from (12.71) that for some c 2 g < oo depending only on D, all x G 3D, z G R™, and all t > 0, 
we have \e tS ^z\ > e~ C29 *|z|. Therefore, for any w G T Xo D, 

> exp ( -c 29 M J WxL^xL *4*4 W \ 



i=0 



c 30 I TTaJ ^ r p ' ' ' 71 x i ^ x i w | • 



(i 



It follows that 

krjl 

n 

i w i 

and, therefore 



1 |7T J ■ ■ -7T J7TJW 

W I T_J • • • 7T_j 7T_j W I 

x i x o 



log |v J I = log I w| + U°g kaJ ■ ■ • 7r 4 7r a J W l ~~ lo S Vx{ ' " 7r 4 7r a: J 1 W l ) 1 

By the Pythagorean theorem, |z| 2 = | z | 2 + (z/|z|, n(x)) 2 |z| 2 . This implies that for some 
C31 < 00, if z G TydD then 

kxz| > (l - c 3 i\x - y\ 2 ) |z|. 
Thus we can find p\ > so small that for some C32 and all \x — y\ < pi and z G T y dD, 

log 1 71-3; z I > log |z| - c 32 \x - y\ 2 . 
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Therefore, 



log |v J '| > log |w| - c 32 \ x l - x fe+il 2l {K-4 +1 |< Pl } ( 3 - 30 ) 



k=l 



7T j ■ ■ ■ 7T j 7T j W | 

X k 1 



+ E 1 (\x j -x 1 ^pi) 10 ^- 



We make p\ smaller, if necessary, so that pi/2 = Ej 2 for some integer j 2 . Note that the set 
of excursions e t j 2 is finite, with cardinality rrij 2 . 

The hitting distribution of 3D for any excursion law H x is absolutely continuous with 
respect to the surface area measure on 3D, because the same is true for Brownian motion. 
This, (13. 4p and Assumption 13.11 imply that with probability 1, for all k = 1,2, ... ,rrij 2 , 
we have |(n(e j 2 (0)), n(e J2 (C — )))l > <5> f° r some random 5 > 0. For large j, because of 
continuity of reflected Brownian motion paths, and because excursions are dense in the 
trajectory, the only points x k+1 such that {x^ — a^. +1 | > p\ can be the endpoints of excursions 

j i = 1; 2, . . . , Ulj 2 . 

i 

Fix a point e^ 2 and let k(j) be such that = e t i 2 - Then x ] k ^_ x — > xj,^ as j — > oo, 
again by the continuity of reflected Brownian motion, and because excursions are dense in 
the trajectory. It follows that for large j, for all pairs (x k ,x k+l ) with \x k — x k+1 \ > p±, we 
have \(n(x k ),n(x k+1 ))\ > 5/2. This implies that, a.s., for some random U > — oo, and all 
sufficiently large j, 

1 7T j ■ ■ ■ 71 J 71 jWl 



E V-v i^} 10 ^^ — ~ - i > u - ( 3 - 31 ) 

At — 1 y X k-1 1 

In view of < K21h and (EOfiT) . for any f3 7 < 1, 



E(E 2 j (3.32) 



< 8E 



mj max sup \xl — X t \\ I + 8E ( |e t j (0) — ((- 

V Y- k - m3 ti+ae{)<t<t{ +1 ,X t £dD J J \k=t 



<c 23£ f + 8E E l e 4( ) - e *i(C 



-/3 7 ^ ^ 

By (13 .4p and (13.111) . the expected number of excursions e s with |e s (C— ) — e s (0)| G [2~ i ~ 1 , 2" 
and s G [0, 1] is bounded by 0332*. It follows that for some C34 < 00, not depending on j, 



E E 1^(0) - e ti^~)\ 2 ^ E c 3i 2- 2l T < c 35 < 00, 

\fc=l / i=ji 

and this combined with (13.321) yields 



supE [Y.Wi -4+x\ 2 ) < °°- 

3>h \i=0 J 
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In view of ( ETSOjl and fl33T]) . 

liminf E(log |v J | —U)> log |w| — limsup E C32 V Ixl — xi, J 2 > —00, 

jwoo \ fe =i / 

so, with probability 1, liminf j_+oo |v J | > 0, and, therefore, |vi| 7^ 0. □ 
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